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Abstract 

The presence of nanoscale inhomogeneities has been experimentally evidenced in several diluted 
magnetic systems, which in turn often leads to interesting physical phenomena. However, a proper 
theoretical understanding of the underlying physics is lacking in most of the cases. Here we 
present a detailed and comprehensive theoretical study of the effects of nanoscale inhomogeneities 
on the temperature dependent spontaneous magnetization in diluted magnetic systems, which 
is found to exhibit an unusual and unconventional behavior. The effects of impurity clustering 
on the magnetization response have hardly been studied until now. We show that nano-sized 
clusters of magnetic impurities can lead to drastic effects on the magnetization compared to that 
of homogeneously diluted compounds. The anomalous nature of the magnetization curves strongly 
depends on the relative concentration of the inhomogeneities as well as the effective range of the 
exchange interactions. In addition we also provide a systematic discussion of the nature of the 
distributions of the local magnetizations. 

PACS numbers: 75.60.Ej, 78.67.Bf, 75.50.Pp 
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Interest in diluted magnetic systems, such as diluted magnetic semiconductors (DMSs)P^ 
and diluted magnetic oxides (DMOs)P^ has continued to surge over the past couple of 
decades owing to their huge potential for spintronics applications. The prime requisite of 
room-temperature ferromagnetism for feasible spintronics devices has led to Curie tempera- 
ture (Tc) being the primary focus of interest in these diluted magnetic materials. Contrary 
to a longstanding belief of defect and inhomogeneity free samples leading to high Curie 
temperatures, recent studies suggest otherwise. Experimental studies have revealed the for- 
mation and existence of magnetic clusters on the nanoscale order, in materials like (Ge,Mn)P 
and (Zn,Co)CP, which in turn gave rise to very high Tc's (>300 K). In fact, very recently we 
have theoretically shown that incorporating magnetic nanoclusters, can lead to drastically 
high critical temperatures (often above room-temperature) in diluted magnetic systems with 
effective short-ranged exchange interactions^. Now, another important aspect which can help 
to further reveal the physical intricacies of these strongly disordered and complex systems 
is the temperature dependence of the spontaneous magnetization M(T). Among the many 
interesting features that magnetization possesses, some worth mentioning are the convex or 
concave nature of the M(T) curves and their critical behavior close to the transition point. 

In the particular case of DMSs, one of the very first observations of ferromagnetism 
in (In,Mn)As^, revealed a surprising non-mean-field like behavior of the spontaneous mag- 
netization with temperature. The experimentally determined magnetization curve had an 
unusual outward concave-like shape which is in stark contrast to the typical convex be- 
havior obtained within the standard Weiss mean-field theory®!, as we \\ as that observed 
in conventional ferromagnetic materials. The (In,Mn)As samples studied in Refj9] were 
reported to be insulating. Similar concave M(T) behavior was also observed in insulating 
samples of Gei_ x Mn x determined by superconducting quantum interference device (SQUID) 
magnetometry and magnetotransport measurements®! Theoretical predictions, based on a 
percolation transition of bound magnetic polarons in the strongly localized regime®', were 
made to explain this non-mean-field like magnetization behavior. Similar magnetization 
behavior in DMS systems, in the insulating regime, was also predicted by other theoretical 
studies based on numerical calculations^^. The deviation in the M(T) behavior from the 
standard Brillouin-function shape was believed to be partly due to the small carrier den- 
sity compared to the localized spin density, as well as due to the wide distribution of the 
exchange interactions and hopping integrals®'. On the other hand, in metallic DMSs, for 
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example Gai_ x .Mn x As for £=0.05—0.10, the magnetization behavior is found to decrease 
almost linearly with temperatur e! 16 * 17 *. This behavior is somewhat intermediate between 
the concave-like M(T) curves in the insulating regime and the classic convex magnetiza- 
tion. However, experimental studies suggest that annealing treatments can have a strong 
effect on the nature of the magnetization behavio r 118 1 19 1 . During annealing the magnetization 
curves are found to become more convex and Brillouin-function-like compared to their lin- 
ear behavior before annealing. This anomalous behavior of the magnetization highlights the 
importance of disorder in these systems. 

However, the effects of correlations in disorder or impurity clustering have hardly been 
considered barring a few cases. In RefJ20l the authors reported that correlated defects lead 
to "mean-field-like" magnetization curves in (Ga,Mn)As. In another theoretical studj^H 
on Co doped ZnO, inhomogeneous phases were shown to be responsible for high Curie 
temperatures in comparison to the homogeneous samples. Now a majority of the existing 
theoretical studies are based on mean-field-like approaches, which is known to be inadequate 
to treat thermal and/or transverse fluctuations and disorder reliably in these systems. In 
RefJ22l the authors have used some complementary theoretical approaches in addition to 
the mean-field theory, which include the dynamical mean-field theory (DMFT), to study 
the temperature dependent magnetization in doped magnetic semiconductors. Nonetheless, 
the effects of clustered defects were not taken into account. The presence of nanoscale 
inhomogeneities can give rise to very interesting and new physics in these diluted systems, 
as was seen in the case of the Curie temperatures^. 

The primary objective of the current manuscript is to investigate the effects of these 
inhomogeneities on the spontaneous magnetization behavior. Here we first study the nature 
of the magnetization in homogeneously diluted systems with no correlation in impurity posi- 
tions, and then extend this to systems containing clusters of magnetic impurities. We observe 
very interesting as well as strong deviations from the homogeneous magnetization curves. A 
non-trivial nature of the spontaneous magnetization in these inhomogeneous diluted systems 
is found to strongly depend on the relative concentration of the inhomogeneities as well as 
the effective range of the magnetic exchange interactions. 

For the sake of simplicity we have assumed here a simple cubic lattice with periodic 
boundary conditions. We fix the total concentration of impurities in the system to x=0.07 
and the total number of impurities in the system is denoted by N imp . The inhomogeneities 
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FIG. 1: (Color online) Illustration of the self-consistent local random phase approximation method 
at a given temperature T. 

are assumed to be of spherical shape of radii tq. The concentration of impurities inside 
each nanosphere is denoted by Xi n . To avoid additional parameters, we restrict ourselves to 
nanospheres of fixed radii r =2a (a is the lattice spacing) and Xj n =0.8. The concentration 
of nanospheres in the system is defined by x ns =Ns/N, where N$ is the total number of 
sites included in all the nanospheres and N=L 3 is the total number of sites. We define a 
variable P/v=(xj n /x)x ns , which denotes the total fraction of impurities contained within the 
nanospheres. For a particular configuration, the nanospheres are distributed in a random 
fashion on the lattice, the only restriction being to avoid any overlap with each other. 

In order to calculate the spontaneous magnetization, we start with the effective diluted 
Heisenberg Hamiltonian describing N imp interacting spins (Sj) randomly distributed on a 
lattice of N sites, given by 

H = - JijPiPjSi ■ Sj (1) 

where the sum ij runs over all sites and the random variable Pi=l if the site is occupied by 
an impurity or otherwise 0. We assume magnetic couplings of the form Jy= Joexp(-|r|/A), 
where r=rj-rj. Both, the choice for such couplings and the values of the damping parameter 
A are discussed below. 

This Hamiltonian (Eq.l) is treated within the self-consistent local random phase ap- 
proximation (SC-LRPA) theory, which is essentially a semi-analytical approach based on 
finite temperature Green's functions. More details on the SC-LRPA theory can be found in 
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Refs J3|23l A schematic illustration of the method is given in Fig. [TJ We briefly summarize 
the method in the following. We start the self-consistent loop with a fully polarized, collinear 
ferromagnetic ground state at T=0 K. The impurity spins are assumed to be classical in 
this case, although the theory is valid for quantum spins as well. Within the SC-LRPA 
formalism, for a given disorder configuration, the local magnetizations at each site (S l z ) 
(i=l,2, ...,Ni mp ) are calculated self-consistently at each temperature. The local magnetiza- 
tion is evaluated using a Callen-like expression^, which relates the local Green's function at 
site % to the local magnetization at this site, 

(S - $0(1 + $Q 2S+1 + (1 + S + $i)^ 25+1 

;i + $i) 2S+1 - $f 

where the local effective magnon occupation number is given by 



/ci\ _ W -^in^Tfj; T^iTUT^;^ 
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We define the retarded Green's function as Gij(u)=f^° oo Gij(t)e luJt dt, where Gij(t)=- 
i9(t)([S^(t), Sj(0)]), which describes the transverse spin fluctuations, and (...) denotes 
the expectation value. By the self-consistent treatment we obtain, for a given tem- 
perature and disorder configuration, the average magnetization which is defined by 
(S™ 9 )= 1 ^Y£S v {S i z)- The accuracy and reliability of the SC-LRPA to treat ther- 
mal/transverse fluctuations and disorder/dilution have been shown on several occasion^ESl. 
Moreover it has also been applied to calculate the Curie temperatures in inhomogeneous 
diluted systems in a recent study^. 

The assumption of the magnetic interactions of the form Jy— Joexp(-|rj-rj|/A) is based 
on the fact that both ab initio studies^^ as well as model calculations^ have shown that 
the exchange couplings in III-V DMSs are relatively short ranged and non-oscillating in 
nature. This non-oscillating behavior results from the existence of a preformed impurity 
band (IB)^. The existence/absence of an IB has been a longstanding debate for several 
years. But in our view, the still often quoted RKKY interactions in these compounds, only 
consistent with the valence band (VB) or perturbative picture, can be ruled out. As shown in 
RefJ25l only the IB scenario could explain the proximity of (Ga,Mn)As to the metal-insulator 
transition as well as the observed red shift in the optical conductivity in this compound. 
Furthermore, recent experimental findings^ definitely corroborate the fact that the Fermi 
level is located within an IB in (Ga,Mn)As. Now, in (Ga,Mn)As, for about 5% Mn, a fit 



of the ab initio exchange couplings provides a value of A of the order of a/2. It should 
be noted that in the case of (Ga,Mn)N the exchange interactions are even shorter ranged. 
In the following we will consider a range of A's, corresponding to relatively long-ranged 
couplings down to shorter ranged ones, and try to analyze their effects on the magnetization 
behavior. In fact the magnetic couplings could also depend on the nanoscale inhomogeneities. 
However, computing these would require extensive calculations involving finite size analyses, 
systematic average over disorder configurations and especially diagonalizing considerably 
large matrices. This is beyond the scope of the present manuscript. Nevertheless, we believe 
that the overall nature of the exchange couplings, especially in wide gap compounds like 
(Ga,Mn)N, will remain essentially unchanged. The interesting feature which we want to 
focus on in the current work is the role the inhomogeneities play in determining the nature 
of the magnetization curves. In the following, the average magnetization is always plotted 
as a function of the reduced temperature T/T*, where T* is the temperature corresponding 
to the case when (S® V9 )=0.001S. We have chosen T* instead of T c as it is difficult to 
determine accurately the critical temperatures from the magnetization curves. However, we 
have checked for some cases, that T* is relatively close to the Tc directly calculated from 
the semi-analytical expression (Eq.(l) in RefJH]). 

To begin with we consider the homogeneously diluted case where the magnetic impurities 
are randomly distributed on the lattice. Fig. [2] shows the average magnetization as a func- 
tion of the reduced temperature T/T*, for different values of A. The magnetization shown 
corresponds to a single configuration of disorder and the system size is L=24. We observe 
that for relatively long-ranged couplings (A=2a), the magnetization curve has a pronounced 
convex shape which is the usual behavior predicted by the mean-field theories of Weiss and 
Stoner^, and observed commonly in conventional ferromagnetic materials. On decreasing A, 
we notice that the convexity decreases and for A=a/3 the magnetization is more linear-like 
over a broad temperature range. In fact similar behavior (linearity) of the magnetization was 
observed in metallic samples of (Ga,Mn)As by magnetic circular dichroism (MCD) studies^. 
This shows that the relatively short ranged interactions are more relevant for the case of 
DMS materials and also vindicates the choice of our exchange couplings. In order to have 
a qualitative idea of the relative change in the behavior of the magnetization with A, we 
have plotted in the inset of Fig. [2] the curvature (k) of the magnetization curves at the spe- 
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FIG. 2: (Color online) Average magnetization as a function of T/T* for different values of A, 
corresponding to the homogeneous case for a single configuration. Inset: Curvature of the {S^"" 9 ) 
curves as a function of A, calculated at T/T*=0.5. The blue dashed line indicates the saturation 
value of the curvature. (Here A=24 3 ). 
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cific value of T/T*=0.5. The curvature is defined by k- 
where u=^. As can be clearly seen, with increase in A the curvature changes significantly, 
increasing by almost a factor of five from A=a/3 to A=2a. For A>2a, one sees that k has 
already saturated and the magnetization has a standard Brillouin shape. This implies that 
A>2a corresponds to the long range coupling regime. 

In Fig. [2j we have shown the magnetization for a single configuration of disorder. However, 
we know that in diluted materials the magnetic properties are often very sensitive to the 
random impurity configurations. Fig. 3 shows the average magnetization calculated for 
25 configurations of disorder corresponding to three different values of A. We have also 
calculated , which gives a measure of the fluctuations of the extremal magnetization 



curves from the configuration averaged one ((Sz Vg )), at the particular value of T/T*=0.6. 



For X=a (Fig. 3(a)), we observe that the spontaneous magnetization is weakly sensitive to 
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FIG. 3: (Color online) Average magnetization for the homogeneous case calculated for 25 config- 
urations, for (a)A=a, (b)A=a/2, and (c)A=a/3. The thick black lines with symbols indicate the 



configuration averaged magnetization, ^=g?|y is a measure of the fluctuation at T/T*=0.6. (Here 
iV=24 3 ). 



the disorder configurations, the overall shape of the curves is unchanged. is found to 

vary within 10% of the configuration averaged magnetization. For X=2a (not shown here), 
these fluctuations were found to be even smaller, varying within less than 5%. For the short 
ranged couplings, X—a/2, the magnetization curves are still convex but the fluctuations are 
stronger now. j^mm is more than doubled as compared to the intermediate range of X—a. 



As can be seen from Fig. 3(b) , some of the curves have a regular convex behavior while some 



are more linear in nature. Now on further reducing A (Fig. 3(c)), the deviations become 
even stronger, and a significant number of configurations exhibit a clear linear temperature 
dependence. The fluctuation at T/T*=0.6 increases by almost a factor four, compared to 
the case of X—a. Even, in some cases, the magnetization profiles are slightly concave toward 
the high temperature. It should be noted that the more linear or concave magnetization 
curves correspond to relatively high Tc's. This figure clearly shows that the disorder effects 
are significantly enhanced in the case of short-ranged interactions. The primary reason is 
that the probability to find regions of weakly interacting impurities increases significantly 
for the case of short-ranged interactions. This is the case in most of the DMSs, where a 
non-trivial behavior of the magnetization is observed. 

So far we have only considered homogeneously diluted systems assuming a fully ran- 
dom distribution of the impurities. Now we move to the case of nanoclusters of magnetic 
impurities. Fig. [4] shows for A=a, the average magnetization of the whole system, (5**°*), 
as a function of T/T*, for four different concentrations of nanospheres, x ns =l%, 3%, 5%, 
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T/T* 



FIG. 4: (Color online) (S^) stands for the total average mag netization ((5*°*)), magnetization inside 
the nanospheres ((S 1 *™)) and outside the nanospheres ((S° ut )). Four different concentrations of 
nanospheres (x ns ): (a) 0.01, (b) 0.03, (c) 0.05, and (d) 0.07, are considered. Pjv is the percentage 
of total impurities contained in the nanospheres. Here A=a, ro=2a, Xj n =0.8 and iV=24 3 . The 
x-axis is in units of T/T* . 

and 7%. In addition, we have also shown the average magnetization inside and outside the 
clusters denoted by (5"* n ) and (S° ut ), respectively. The curves shown here correspond to a 
single configuration of disorder, the variation with disorder configurations will be discussed 
in what follows. We immediately observe that in the presence of inhomogeneities the spon- 
taneous magnetization has a non-trivial behavior and exhibits a drastically different nature 
when compared to the homogeneous case (Fig. [2]). This can be clearly seen even for the 
lowest concentration of nanospheres. For x ns =0.01, for which 11% of the total impurities 
are inside the nanospheres, (Sl ot ) decreases rapidly till about T/T*~0.5, then it becomes 
concave and decays slowly toward the higher temperatures. By gradually increasing the con- 
centration of the nanospheres, an interesting change in the average magnetization behavior 
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is observed. For x ns =3% (P/v=0.34), (Si 01 ) falls off less sharply at low temperature, for 5% 
it is almost linear over the entire temperature range, and for 7% it becomes more convex. 
Thus a crossover in the curvature of (5"*°*) appears at x„ s ~0.05. On the other hand, (S l z n ) 
exhibits a clear convex nature which does not change with x ns . This indicates that the av- 
erage magnetization inside the clusters remains almost uniform and is mainly controlled by 
the intra-cluster couplings. We can clearly see that the inhomogeneities have a very strong 
effect on the impurities outside the clusters. (S° ut ) has a very pronounced concave nature 
which can even be seen for relatively small x ns . The slope at low temperatures becomes 
steeper with increasing concentration of nanospheres. For example, at T/T*~0.3, (S° ut ) 
has a value of 0.85 for x ns =0.01, 0.62 for x ns =0.03, and about 0.4 for x ns =0.05. Similar 
concave behavior of the temperature dependent magnetization is observed in the case of 
some insulating DMS materials^. However, in most of the cases studied until now clustering 
effects have hardly been considered. Note that we have also performed the calculations for 
A=2a (not shown here). In this extended coupling regime, it was found that the effect of 
inhomogeneities are very weak. (5"*°'), (S l z n ), and (S° ut ) exhibit a convex nature and were 
found to be relatively close to each other. 

In an earlier study, we have shown that relatively short-ranged couplings appear to have 
spectacular effects on the Curie temperatures in the presence of nanoscale inhomogeneities^. 
Short-ranged interactions are more relevant from the practical point of view. As men- 
tioned above, the magnetic couplings in most III-V DMS materials ((Ga,Mn)As, (Ga,Mn)N, 
(Ga,Mn)P,...) are effectively short-ranged in nature. In Fig. [5] we show (5*°*), (S l z n ), and 
(S° ut ) as a function of T/T*, for four different concentrations of nanospheres, in the case of 
relatively short-ranged couplings, namely X=a/2. To start with, we first discuss the results 
for a single configuration of disorder. For the lowest x ns (Fig. [5Fa)), the behavior of (Sl ot ) 
is almost similar to that observed in the case of X=a (Fig.gja)). But on increasing x ns fur- 
ther (Fig. |5^b)), we immediately observe that (Si 01 ) decreases much more rapidly at lower 
temperatures followed by a slow decay toward the high temperatures. Also, in Fig. |5](c) 
and (d), an inflection appears in (Sl ot ) around T/T*~0.6. In this case (S z n ) too exhibits a 
non-trivial behavior for all values of x ns , which is unlike the case of X=a. For example, for 
x ns =0.05 (Fig. |5^c)), there is a shoulder-like feature in (S l z n ) around T/T*~0.05, which is 
absent for \=a (Fig. |4|c)). Thus unlike the case of X=a, where the intra-cluster couplings 
dominate, there are other relevant couplings, like the inter-cluster ones and those between 
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FIG. 5: (Color online) (S^) stands for the total average magnetization ((S"* ')), magnetization 
inside the nanospheres ({S l z n }) and outside the nanospheres ((S° ut )). Four different concentrations 
of nanospheres (x ns ): (a) 0.01, (b) 0.03, (c) 0.05, and (d) 0.07, are considered. Pjy is the percentage 
of total impurities contained in the nanospheres. Here A=a/2, ro=2o, X{ n = 0.8 and iV=24 3 . The 
x-axis is in units of T/T* . 

the cluster impurities and bulk impurities, which come into play. On the other hand, the 
{S° ut ) curves are typically concave for all considered x ns , and exhibit a long tail toward the 
higher temperatures. They exhibit a sharp fall-off at low temperatures with increasing x ns . 
At T/T*~0.2, for x ns =0.01, the value of (S° ut ) is about 0.8 which falls rapidly to almost 
0.3 for x ns =0.05. With increasing x ns , the concentration of impurities outside the clusters 
gradually decreases, leading to an increase of the typical distance between them. Conse- 
quently the impurities outside interact more weakly with each other and this explains the 
sharp fall-off in (S° ut ) at lower temperatures. 

In the previous two figures we have discussed the cases for a single configuration of disor- 
der only. Let us now analyze how the results depend on the random cluster configurations. 
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FIG. 6: (Color online) Total average magnetization as a function of T/T* for x ns =0.05 calculated 
for 50 configurations (thin continuous lines), for (a)A=a, and (b)A=a/2. The thick black lines 
with symbols indicate the configuration averaged magnetization (S 1 * *). ^p|y is a measure of the 
fluctuation at T/T* =0.6. (Here r =2a, x in =0.8 and iV=24 3 ). 

Fig. 6 shows the average magnetization calculated for 50 configurations of disorder, for 
x ns =0.05. We compare the case of the intermediate couplings to the short-ranged ones. One 



sees for the case of X=a (Fig. 6(a)), that the magnetization curves have an almost linear 
similar shape over the entire temperature range. In fact the configuration averaged magneti- 
zation can be well approximated by (5'* *)~(1 —T/T*). On the other hand, the configuration 



averaged magnetization has a pronounced concave nature for X=a/2 (Fig. 6(b)). The de- 
cay slope at low temperatures (T/T*<0.2) is twice that of X=a. Concerning the disorder 
fluctuations, one clearly sees that they are much stronger in the case of the short-ranged 
couplings. For example, the fluctuation of (S 1 " 1 ) at T/T*=0.6 is found to be more than 
doubled compared to that of X=a. For X=a/2, most of the magnetization curves are con- 
cave in nature, while some exhibit a linear-like behavior. The typical separations between 
the clusters play a decisive role, in the presence of short-ranged couplings. This has been 
discussed in the context of Curie temperatures in inhomogeneous systems^. Here it should 
be noted that the concave-like curves correspond to higher critical temperatures, while the 
linear ones coincide with relatively low Xc's. Finally, we have found that for relatively ex- 
tended couplings (A~2a) (not shown here), (S 1 " 1 ) exhibits a more convex-like behavior. The 
fluctuations resulting from the disorder configurations were found to be much smaller than 
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FIG. 7: (Color online) Configuration averaged magnetization (Sl ot ) for different system sizes N=L 3 , 
corresponding to (a) X=a, and (c) X=a/2. The total number of configurations is of the order of 
few hundred. The respective fluctuations 7==- at T/T*=0.6, as a function of 1/N is shown in (b) 
for X=a, and (d) for A=a/2. (Here x ns =0.05, r$=2a, and Xi n =0.8). 

that of X=a. 

We now proceed further with the study of the finite size effects. Indeed, in inhomogeneous 
systems one naturally expects finite size effects to be much stronger than in homogeneously 
diluted ones. These effects can be even more drastic for larger size of inhomogeneities. The 
configuration averaged total magnetizations as a function of T/T*, corresponding to X=a 
and X=a/2, are shown respectively in Fig. |7[a) and (c). The calculations are performed 
over different system sizes varying from N=12 3 up to iV=28 3 . (S 1 * *) is the averaged mag- 
netization obtained over a sample of few hundred disorder configurations. In both cases we 
observe that the magnetization curves are very similar and the finite size effects are very 
weak. Let us now discuss the size dependence of the fluctuations, j=- , (with respect to 
disorder configurations) of the average magnetization (Si 01 ). Fig. mb) and (d) show j=- 
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FIG. 8: (Color online) Distribution of local magnetizations (5*) for A=a/2, at three different 
temperatures T\, T2, and T3, which correspond to (5*°*) =0.755', 0.45 and 0.15 respectively. Left 
column: Inhomogeneous case for x ns =0.05 (with tq=2cl and Xj n =0.8). Right column: Homogeneous 
case. (The red dashed lines indicate the values of (5*°')/5 corresponding to the three different 
temperatures). 

at T/T*=0.6 as a function of 1/N, corresponding to X=a and X—a/2, respectively. For 
X=a, it is almost constant by varying the system sizes, with a value of around 0.42. In- 
terestingly for the short-ranged couplings, j=- decays monotonously, and indicates that 

Wz° ) 

it saturates to a value of about 0.7 in the thermodynamic limit. It is important to note 
that for homogeneously diluted systems t=t vanishes in the thermodynamic limit. This 

(S z ° ) 

means that the total average magnetization is not a self averaging quantity in the presence 
of inhomogeneities. 

In order to get a deeper insight into the temperature dependence of the magnetization 
in the presence of inhomogeneities, we analyze the nature of the distributions of the local 
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magnetizations. In Fig. |8j we compare the distributions of (S^'s in the inhomogeneous 
systems with those of the homogeneously diluted ones, for X—a/2. The distributions are 
shown for three different temperatures 7\, T 2 , and T 3 , which correspond to (S t z ot )=0.75S , 
OAS, and 0.1S, respectively. The distributions at the lowest temperature 7\ are almost 
similar for both cases, except for an extended tail toward the small magnetization values 
in the inhomogeneous system. At the intermediate temperature T 2 , the distributions reveal 
a more pronounced difference. It has a Gaussian-like shape for the homogeneous systems, 
whereas it has a clear bimodal nature in the inhomogeneous ones. It is also wider than 
that of the homogeneous case. There is a significantly higher weight below (5**) =0.2 for 
the inhomogeneous systems. This relatively high weight corresponds to the impurities out- 
side and far from the clusters. The broad peak around (S"*)~0.7 is associated with the 
impurities inside the clusters. Finally, for the highest temperature T 3 , the distribution in 
the homogeneous systems is unimodal with a small tail at large (-S^'s and the half-width 
is of the order of (S'* o *)=0.1. On the other hand, in the inhomogeneous compounds, we 
observe a very strong weight for (S'*)<0.05 and an extended tail which goes up to 4(5"*°*). 
For T=T3, an estimate of the percentage of impurities within the range of 0.05<(5'*)<0.15 
gives 65% for the homogeneous systems and only 15% in the presence of inhomogeneities. 
To conclude this discussion, we provide in Fig. 9 a real space illustration of the variation 
of the local magnetizations. The 2D-snapshots, are shown for both kinds of systems, at the 
temperatures T±, T 2 , and T 3 , corresponding to X=a/2. 

Thus, to summarize, we have presented a detailed and comprehensive study of the effects 
of nanoscale inhomogeneities on the spontaneous magnetization in diluted magnetic systems, 
which had hitherto been unexplored. We have compared the average magnetization behavior 
in inhomogeneous systems to that of the homogeneously diluted case, for different ranges 
of the magnetic couplings. Unlike the convex nature of the spontaneous magnetization in 
homogeneous systems a linear temperature dependence (over the entire temperature range) 
is obtained in the inhomogeneous compounds for intermediate couplings. Whereas it ex- 
hibits a pronounced concave shape in systems with short-ranged interactions. Additionally, 
the local magnetizations show bimodal and wider distributions in contrast to that of the 
homogeneously diluted systems, where it remains unimodal. The finite size analyses have 
revealed that the fluctuations of the average magnetization (with respect to disorder con- 
figurations) remain finite in the thermodynamic limit in inhomogeneous systems, unlike the 
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FIG. 9: Snapshots of the local magnetizations in a 2D-plane for X=a/2 . The top panel corresponds 
to the homogeneous case, and the bottom panel to the inhomogeneous case for x ns =0.05. Ti, T2, 
and T3 denote the temperatures when (S t z ot )=0.75S, OAS and 0.15 respectively. (Here iV=24 3 ). 

homogeneously diluted case for which it vanishes. We believe it would be of great interest 
to corroborate our findings by future experimental studies. 
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